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We investigate the temporal evolution and spatial propagation of branching annihilating 
random walks in one dimension. Depending on the branching and annihilation rates, a 
few-particle initial state can evolve to a propagating finite density wave, or extinction may 
occur, in which the number of particles vanishes in the long-time limit. The number parity 
conserving case where 2-offspring are produced in each branching event can be solved ex- 
actly for unit reaction probability, from which qualitative features of the transition between 
propagation and extinction, as well as intriguing parity-specific effects are elucidated. An 
approximate analysis is developed to treat this transition for general BAW processes. A 
scaling description suggests that the critical exponents which describe the vanishing of the 
particle density at the transition are unrelated to those of conventional models, such as 
Q\ . Reggeon Field Theory. 
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1: Introduction 

> ■ 

In the branching annihilating random walk (BAW), a single random walk branches at some 
specified rate and two random walkers annihilate at another rate when they meet. 1-4 As 
a function of these rates, the number of random walkers may grow without bound, reach 
a finite limiting number, or vanish asymptotically. Our goal, in this paper, is to determine 
some of the long-time properties of this BAW process. We are particularly interested in 
understanding the kinetics and density distribution when the initial state consists of a 
small number of localized particles. 

Interest in this process has several motivations. First, considerable theoretical ef- 
fort has been devoted to establishing the existence of and quantifying the non-equilibrium 
phase transition between "propagation" and "extinction" for a variety of interacting par- 
ticle systems 5 which are closely related to BAWs. Here the term extinction refers to the 
situation where annihilation dominates over branching and an initially localized popula- 
tion of particles ultimately disappears. In the complementary case, branching dominates 
over annihilation and an initially localized population evolves into a propagating wavefront 
which advances into the otherwise empty system. Typical examples of these phenomena 
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include the contact process, 6 Schlogl models, 7 as well as directed percolation and Reggeon 
field theory. 8 The propagation phenomenon is also a discrete realization of the "Fisher 
wave", 9 ' 10 which describes the continuum dynamics of an initially localized single-species 
population whose evolution is influenced by diffusion, as well as by both (linear) birth 
and (quadratic) death mechanisms. The relation between the continuum description of 
the Fisher wave and the corresponding BAW is tenuous and comprehensive investigations 
of discrete BAW models would be helpful to understand better the relation with the con- 
tinuum counterpart. Second, there is a direct correspondence between the 2-offspring 
BAW and the interface dynamics in the reaction-limited monomer-monomer surface re- 
action model in one dimension. 11 ' 12 For the surface reaction, a lattice is filled with A 
and B particles and the ensuing dynamics is defined by randomly and repeatedly selecting 
a nearest-neighbor AB pair and changing it to either AA, BB, AB, or BA at specified 
rates. The dynamics of AB interfaces is identical to that of the individual particles in the 
2-offspring BAW. 

In addition to connections with various non-equilibrium systems, the BAW model 
is amenable to theoretical analysis. Somewhat surprisingly, we find that the transition 
between propagation and extinction is controlled by detailed features of the underlying 
discrete process. In particular, the exact solution of the 2-offspring BAW model in one 
dimension reveals that propagation occurs only for infinite branching rate and the parity 
of the initial number of particles essentially influences the long-time kinetics. 

In the next section, we outline several basic facts about the BAW process. The general 
conditions which lead to propagation or extinction are discussed. In Section 3, we present 
an exact solution for the evolution of the 2-offspring BAW process A — > 3A and 2A — > 0, 
in the case where the probability of reaction when two particles meet, k, is unity. For this 
case, extinction is paradoxically found to occur for all non-zero values of the diffusion rate. 
When there is only branching, a behavior intermediate to propagation and extinction 
occurs. In Section 4, we present an approximate description for the transition between 
propagation and extinction by solving a truncated hierarchy of rate equations for multi- 
particle correlation functions. The primary result of this treatment is that the transition 
emerges naturally at the next level of approximation beyond mean-field theory. In Section 
5, we present numerical simulation results to support our various theoretical predictions, 
and then conclude in Section 6. 



2: Models and Phenomenology 

We define the branching annihilating random walk (BAW) on a lattice as follows. A particle 
is picked at random. It can undergo either nearest-neighbor diffusion or it may branch 
with respective rates D and r. In a diffusion attempt, a random direction is picked and 
the particle moves to the target site if it is unoccupied. If the target site is occupied, then 
annihilation of the incident and target particles occurs with probability k. Otherwise, 
the incident particle remains at its original position. The details of the branching step 
depend on the number of offspring produced. In the 2-offspring BAW, if a branching 
attempt is selected (at rate r), then branching to the two nearest-neighbor sites occurs 
with probability 1 if both neighbors are empty. If one or both neighbors are occupied, 
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then branching to both neighbors occurs with probability k. In this branching, if a newly 
created particle is placed on a previously occupied site, then both particles are removed. 
An analogous procedure is employed for the 1-offspring process. Namely, in a branching 
event, if the target is occupied, then branching occurs with probability k. This branching 
entails immediate annihilation of the newly-created and target particles. While these 
microscopic rules are somewhat involved when both branching and a finite reaction rate 
are operative, any reasonable discrete realization of the continuum process is anticipated 
to give qualitatively similar results. (An exception is the case of parallel dynamics, as 
opposed to the serial updating assumed here.) Since D and r have dimensions of inverse 
time and k is dimensionless, two dimensionless parameters which characterize the system 
are D/r and k. 

As mentioned above, the BAW process with an initially localized particle population 
generically exhibits a transition from extinction (which is of one of two types; see below) 
to propagation. However, this transition appears to be at odds with a standard mean-field 
treatment of the model. Indeed, such a description of the BAW leads to the reaction- 
diffusion equation 

= ap(r, t) - bp(r, t) 2 + DAp(r, t), (2.1) 

where p(f, t) is the particle density at position r at time t. Here a, 6, and D correspond 
roughly to the branching rate, the reaction rate, and the diffusion rate, respectively. This 
equation admits a propagating Fisher wave solution of asymptotic density p^ = a/b which 
advances into the low-density region with a velocity proportional to y/aD, whenever a > 0. 
Since a is implicitly positive in this description, there appears to be no mechanism for 
extinction. 

To generate a transition to extinction, a mechanism which changes the sign of a is 
needed. Such a mechanism, however, is easily realizable. At low densities, mean-field 
theory erroneously postulates that the encounter probability varies as the square of the 
particle density, since all spatial correlations are neglected. However, a more complete 
treatment of correlations leads to an encounter probability which is linear in the particle 
density in the low density limit. This follow from the exact reaction-diffusion equation 

do(f t) 

J ^- L = aexactp(r, t) - 6exactP2(0, t) + DAp(f, t), (2.2) 

where P2(0, t) denotes the two-point correlation function at zero separation, and branching 
is accounted for by a exa ct > 0. However, for p <C \J a exact /D, any reasonable theory must 
give p2(0, t) ~ (p) + (p) 2 because the birth mechanism produces nearby pairs of particles. 
Thus, extinction may arise in low spatial dimension because branching and annihilation 
terms are of comparable magnitude in Eq. (2.1); in particular, the sign of the coarse-grained 
branching rate a = a ex act — ^exact may change as external parameters are varied. In the 
following, we shall use r for the branching rate and k for the annihilation probability. 

In numerical simulations, a transition from propagation to extinction generically oc- 
curs as the ratio of diffusion to branching increases (Fig. 1). As the annihilation probability 
(equivalent to the reaction rate) increases, the extinction transition is shifted to lower val- 
ues of D/r. A crucial determinant for extinction is the "parity" of the branching event. 
For example, an even-offspring BAW process which starts with an odd number of particles 
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cannot exhibit true extinction. Rather, the number of particles remains bounded when 
D/r is large. Thus extinction comes in two versions: In non-parity preserving models 
(odd-offspring BAW), or in parity preserving models with an even number of particles 
initially, there is true extinction, whereas in the complementary cases, the particle number 
remains bounded. 

The special case of the 2-offspring BAW in the limits of pure branching, D/r = 0, 
and also complete reaction, k = 1, deserves emphasis because of its peculiar and easily 
visualized features (Fig. 2). The crucial aspect of this situation is that a pair of nearest- 
neighbor particles diffuses rigidly under the action of the branching process. Through the 
branching process, a single initial particle spawns rigid pairs at a finite rate which then 
diffuse freely. In this situation, the number of particles grows as \ft. It is only when there 
is non-zero diffusion that extinction can occur. 



3: Exact Solution for the 2-Offspring BAW in One Dimension 

In this section, we derive the exact solution for the one-dimensional 2-offspring BAW 
in the case of unit reaction probability. This provides a complete description of the tran- 
sition between extinction and propagation. Our approach is based on solving for the time 
evolution of Pj^(t), the probability that there is an odd number of particles in the interval 
[j, k}. (Notice that the particle density at site j is simply Pjj(t).) In a spatially homo- 
geneous system, Pj^(t) is a function of j — k only and satisfies a soluble single- variable 
diffusion equation 12 while for the heterogeneous system, Pj,k(t) is a function of both j and 
k. In the continuum limit, the master equation for Pj ; k(t) can De written as a 2-dimensional 
diffusion equation with a radiation boundary condition. 

To write the master equation for Pj^t), define Dq and ro as the hopping and branch- 
ing rates, respectively, in a elemental event. In a single time step, the parity of the number 
of particles within [j, k] can change only by hopping or branching events at j — 1, j, k, or 
k + 1. From straightforward but somewhat tedious bookkeeping of all these microscopic 
processes (Fig. 3), we obtain the following equations for P J)fc for j > k + 1, 

Pj, k (t+ 1/N) -Pj, k (t) = £°±I°(P j+ltk ( t ) + Pj, k+1 (t) + P jt fc _i(f) + P j -i,k(t) -4Pj, k(t)), 

(3.1) 

while for j = k, 

Pjj(t + 1/N) - P hJ {t) = £°±I°( Pj+lij (t) + Pj, j+1 (t) - 2P iJ (f)) - ^Py(f). (3.2) 
The first equation can be extended to include j = k by introducing the boundary condition, 

(Do + roXPfj-iCf) + P j+ ij(*)) = *roPjj(t). (3.3) 

Here we have taken the time increment for an individual birth or diffusion event to be of 
the order of the inverse system size so that each particle will typically be updated once in 
a unit time step. 
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These equations can be simplified by transforming to the coordinates x = j + k and 
y = j — k. After taking continuous time derivatives, one finds, for y > 1, 

Px, v (t) = (^o+r )(Px+i, y +i(t)+Px+i, y -i(t)+Px-i, y +i(t)+Px-i, y -i(t)-4P :!; , y (t)), (3.4) 
while for y = one has, 

Px,o{t) = ( D o + r )(P*+i,i(*) + P x -i,i{t) ~ 2Px,oW) " 2D P xfi (t). (3.5) 
These two equations can again be unified by imposing the boundary condition, 

(D + r )(P x+1 ,. 1 (t) + P x -i,-i{t) ~ 2 ^,o(t)) = -2D P Xj0 (t). (3.6) 

To obtain a non-singular continuum limit for these master equations, (Do + ro) must 
be of order (Ax)~ 2 , whereas from the boundary condition, (Do + ro)/ Do should be of order 
(Ax) -1 . These two restrictions, in turn, imply that r ~ (Ax) -2 and D ~ (Ax) -1 . That 
is, an infinitely large microscopic birth rate is needed to counteract the instant recom- 
bination of newly- formed pairs in the continuum limit, thus ensuring a finite continuum 
birth rate. In the following, we simply replace the microscopic combination Dq + ro by 
ro — since ro is infinitely larger than Do — and then take the continuum limit. Denoting the 
continuum limits of (Ax)D and (Ax) 2 ro as D and r, respectively, the equation of motion 
for the interval occupation probability P(x,y) becomes 

dP(x,y) 2 

— _ — =rS7 P(x,y), (3.7) 

with the corresponding boundary condition 

dP(x,y) D 

— — \ V =Q = —P{x,y)\ v=0 . (3.8) 

When D/r = 0, a slightly different line of reasoning must be used to determine the equation 
of motion, but the diffusion equation given above remains valid. 

This equation can be solved by writing the solution as a superposition of incoming and 
outgoing modes in y in Fourier space, and imposing the boundary condition to relate the 
amplitudes of the incoming and outgoing modes. Upon inverting the Fourier transform, 
one finds 

P(x^xo^t) = [ ^I^e'M^oJ-rOfcS+fcS)* ( e ik 2 (y-y ) + * 2 + iD/r e - ik2(y+yo) \ 
J (27r) 2 V k 2 -iD/r J 

(3.9) 

where P(x, y; xo, yo', t) denotes the propagator from (xo, yo) to (x, y) in time t. For comput- 
ing the total number of particles, we only need the propagator with final spatial co-ordinate 
y = 0. In the long time limit, this restricted propagator has the form 

nxAx ,y ;t) * -£-^ eM Jj-^l±A h (3.10) 



5 



which arises by neglecting the dimensionless "radiation" length r/D compared to the 
diffusion length VDt. 

For a single particle present initially at the origin, P(xo, yo; t = 0) is zero for y < \xo\ 
and is equal to one for yo > \xq\. This corresponds to P(xo,yo;t = 0) = Q(yo — \xq\), 
where Q(x) is the Heaviside step function. Thus for the density distribution we find 

poo j-yo 

c(x,t)= / dy / dx o P(x,0;x o ,yot), 
Jo J -yo 

e -x 2 /(8rt) ^3-11) 



(D/r)V2nri 



For the case of two particles which are originally separated by a distance d, P(xo,yo',t = 
0) = 1 when | \xq\ — yo\ < d, and P(xo,yo;t = 0) = otherwise. This initial condition 
leads to an expression of a similar form to Eq. (3.11), except that there is an additional 
multiplicative factor of d/y/ri. Thus the total number of particles, N(t) = J c(x,t) dx — > 
2r/D as t — > oo when one particle is initially present, and N(t) ~ (rt) -1 / 2 when two 
particles are initially present. Generally, an initially localized group with an odd number of 
particles leads to a finite number of particles as t — > oo, while a localized initial population 
with an even number of particles leads to N(t) oc t -1 / 2 . 

If there is no diffusion, then Eq. (3.8) reduces to a Neumann boundary condition. 
With this simplification, one finds the following density distribution for a single particle 
initial condition, after several straightforward steps and without any approximations other 
than those involved in the continuum limit, 



poo py 

c(x,t)= / dy / dx P(x,0;xo,yo;t), 

J0 J-yo 



( 3 - 12 ) 

1 / - -2 I x 



1 - erf 



4 V \2^/rt 



Thus for a single particle initial condition, N(t) oc intermediate to the limiting cases 
of Fisher wave propagation, where the particle number grows linearly in t, and extinction, 
where the particle number either remains constant (through parity effects) or else decays 
as t -1 / 2 . Another interesting feature of the distribution is that it qualitatively resembles a 
Gaussian. This can be made plausible by consideration of the underlying discrete process. 
When there is only branching, it is easy to verify that a nearest-neighbor particle pair 
propagates diffusively as a rigid unit (Fig. 2). When two such soliton-like excitations 
meet, they merely "scatter" without any change in their form. As a single initial particle 
evolves, there is production of pairs at a finite rate which then diffuse freely within a region 
of length \fri. Conversely for an even number of particles initially, "interference" between 
the offspring of the initial seed particles leads to the density vanishing as t -1 / 2 in the 
long-time limit. 

When the reaction rate k < 1, then the transition between extinction and Fisher wave 
propagation can be observed numerically, and the transition line between these two states 
has the qualitative form in the D/r - k phase plane indicated in Fig. 1. Our exact analysis 
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also provides a basis for a scaling description of this transition. For the total number of 
particles, we make the following scaling ansatz, 



N(t) ~ t^<&(et a ), (3.13) 

where e(k) is the deviation of D/r from its critical value. This parameter will be considered 
positive in the region where propagation occurs and negative otherwise. From the fact that 
N(t) grows linearly with time for e > and decays as t -1 / 2 otherwise (for the more generic 
case of parity non-conserving dynamics or an even number of initial particles and parity 
conserving dynamics), one obtains 

~ x (1 -^ /a (x>0), 

(3.14) 

„ | a .|-(l+2/9)/(2a) {x<Q y ^ > 



Since the exact solution shows that N(t) remains bounded as t — > oo for e = 0, (3 must 
equal 0. (As above, we consider the case where parity is not conserved, or that the initial 
number of particles is even for parity conserving dynamics). Further, for small values of 
D/r, this ratio appears in the exact solution in combinations of the form Dy/i/r, implying 
that a = 1/2. It therefore follows that N(t) oc e 2 t for small positive e. This cannot be 
checked against our exact solution, since the e > regime cannot be reached when k = 1. 
On the other hand, for negative e, the behavior N(t) oc ( |e| v^) — 1 is predicted, in agreement 
with Eq. (3.11). 



4: Factorization of the Multiparticle Rate Equations 



As mentioned in Section 2, the single-particle reaction-diffusion equation cannot ac- 
count for the transition between extinction and propagation unless there is a sign change 
in the coefficient of the linear term in the concentration. Although such a sign change 
can be justified heuristically, it is worthwhile to present a systematic continuum approach 
which leads to this mechanism. Our approximate treatment of the multiparticle rate equa- 
tions accomplishes this task. This approach also has the advantage that it can be applied 
straightforwardly to different microscopic branching and reaction mechanisms. In contrast, 
while the exact solution of the last section provides a complete description of the transition 
in a special case, this method is neither generalizable nor physically intuitive. Given the 
nature of the approximation in the multiparticle rate equations, we anticipate that our 
results will not depend quantitatively on the number of offspring in a branching event 
(except for parity-specific features), but rather, will be generic to the transition between 
extinction and propagation. 

For simplicity, we study the 1-offspring BAW in one dimension which can be repre- 
sented as 

A^2A A + A^O. (4.1) 

It should be mentioned that the 1-offspring BAW with a unit reaction probability has been 
investigated by rigorous mathematical techniques. 1 These approaches have established the 
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existence of a transition between propagation and extinction and provided weak bounds 
for the critical value of D/r. Our approximate method locates the transition for all val- 
ues of the reaction probability. We first determine the first two in the hierarchy of rate 
equations for the multiparticle densities in this process. These equations will be closed 
by factorizing three-particle densities as products of two-particle densities. To write the 
hierarchy of rate equations, define p^j as the probability that the set {k} is occupied. 
Thus po is the probability that site is occupied, po,i is the probability that sites and 
i are simultaneously occupied, Po,ij is the probability that 0, i, and j are simultaneously 
occupied, etc.. The rate equation for po is found by enumerating all configurations in 
which an elemental event changes the occupancy of site 0. For example, if site is empty 
and site 1 is occupied (which occurs with probability po — po,i); then by either branching 
or hopping to the left, site can become occupied. There is a similar contribution if site 
— 1 is occupied. Similarly, there are three elemental events which lead to a decrease in the 
occupancy of site 0, as illustrated in Fig. 4(a). Summing these contributions, leads to the 
rate equation 

Po = rp -[r+ (2D + r)k]p 0tl (4.2a) 

Similarly, the rate equation for po,i is obtained by enumerating all 3-site configurations 
for which an elemental event changes the occupancy of sites or 1. These configurations 
and the associated rates of the processes which change po,i are shown in Fig. 4(b) and lead 
to 

Po,i = rp - (D + r) (1 + fc)p ,i + (D + r)p , 2 - [r + (2D + r)k\p w (4.26) 
The equations for p ,i for z > 2 are obtained similarly (Fig. 4(c)), 

p 0}i = (D + r)p ,i-i - 2-Dp ,i + (D + r)p ,i+i ~[r+ (2D + r)fc](po,t,i+i + Po,i-i,i)- (4.2c) 

To be soluble, these exact equations must now be closed by a suitable truncation. A 
standard approach is to truncate the equations at the 2-particle level by the ansatz 

Po,i,i+i = Po,iPi,i+i/po = Po,iPo,i/Po- (4-3) 

While heuristic, this approximation yields reasonable results and turns out to be simpler 
than the standard Kirkwood truncation. 

To establish the existence of propagation, we investigate the existence of a non-zero 
time-independent solution to the truncated rate equations. The transition between propa- 
gation and extinction is then identified by the locus where this steady-state density vanishes 
as k and D/r are varied. To streamline the notation, let x = D/r and p(i) = po,i/Po- The 
steady-state rate equations reduce to 

l-[l + (l + 2x)fc]p(l) = 0, 

1 - (1 + x)(l + fc)p(l) + (1 + x)p{2) - [1 + (1 + 2x)k]p(l) 2 = 0, 

(1 + x)p(i - 1) - 2xp(i) + (1 + x)p(i + 1) - [1 + (1 + 2x)k]p(l)(p(i) + p(i - 1)) = 0. 

(4.4) 

From the first of these equations, 



Using this in the second equation then gives 



P(2) = (l + fc + 2xfc)(l + ,) - (46) 

Substituting the static value for p(l) into the last of the rate equations, a constant coef- 
ficient recursion equation for p(i) results. Using the boundary conditions supplied by the 
equations for p(l) and p(2), the solution is, 

sck { x "\ ^ 

" w = " + tt*T2^(tt^J • (47) 

with 

1 — x/c , 

"=mT2S- (48) 

Here p = lim^oo p(i) is just the equilibrium single-particle density. Thus a positive 
solution for p exists only when x < 1/k, corresponding to the propagating phase of the 
1-offspring BAW. When x = x c = 1/k, the equilibrium density vanishes; this defines the 
phase boundary between propagation and extinction. Note further that the decay of the 
equilibrium concentration is linear in e as e = x — x c — > + . If the scaling predictions of 
the preceding section are correct, the particle number should grow as e 2 t, for e — > 0~. On 
the other hand, the total particle number is the particle density multiplied by the width 
of the propagating wave. This latter quantity should grow as et, that is, the front velocity 
goes linearly to zero as the critical value of x is approached. 



5: Numerical Simulations 



To confirm the above analytical results and visualize the evolution of the system, we 
have performed numerical simulations of the 1- and 2-offspring BAW in one dimension 
according to the rules outlined in Section 2. For concreteness, we have fixed the hopping 
rate to be unity so that the variables in the simulation are the branching rate r and the 
reaction probability k. 

First consider the 1-offspring BAW for which it is known 1 ' 2 that a transition between 
propagation and extinction occurs for a non-zero value of D/r. Starting with a single 
particle at the origin and with the parameterizations of our lattice model, the transition 
occurs at x = D/r = 1/.89. For slightly larger values of x, the number of particles 
initially increases, but ultimately decays to zero. For example, for x = 1/.88, the average 
number of particles gradually increases to 4.8 for t ~ 400, but then decreases to for 
longer times. By t = 1.5 23 ~ 11223, only about .1% of the initial configurations are still 
active. For 400 < t < 10, 000, the spatial distribution of the ensemble of surviving particles 
appears visually to be well approximated by a Gaussian (Fig. 5(a)). In this time range, 
the reduced moments of the spatial distribution are = (x 4 ) / (x 2 ) 2 = 2.8 - 3.02 and 
1716 = (x 6 )/(x 2 ) 3 = 13.0 - 13.9, while for a Gaussian, the corresponding values are = 3 
and m 6 = 15. This behavior suggests that the effect of branching is irrelevant in the scaling 
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sense, and that the spatial evolution of the 1-offspring BAW coincides with that of a single 
purely random walk. On the other hand for x = 1/.95, the simulations clearly show that 
a Fisher-like wavefront forms which then propagates at a finite velocity (Fig. 5(b)). 

For the 2-offspring BAW, we have approximately mapped out the phase boundary 
between extinction and propagation (Fig. 1(b)). Simulations clearly indicate that the 
phase boundary intersects the line k = la,tD/r = 0. Thus for k = 1 there is extinction 
for all D/r except when D/r = 0. At this special point, our exact solution showed that 
the number of particles grows as \ft for a single particle initial state. Another amusing 
feature of the 2-offspring BAW that we have verified numerically is the essential dependence 
on parity. For example, if the initial state consists of two widely separated particles, then 
propagation associated with two independent particles is observed for early times, followed 
ultimately by a decay in which the number of particle goes to zero. 

Finally, in the course of simulations we discovered an interesting difference between 
serial and parallel updating. When all lattice sites are updated simultaneously, the 2- 
offspring BAW exhibits propagation, that is, a Fisher wave, even for finite values of D/r. 
Thus, the BAW process is an example of a system where the method of updating, whether 
parallel or serial, leads to non-trivial differences in the resulting dynamics. We chose to 
focus on serial updating, since it is more physically motivated. 

6: Conclusions 

In summary, we have analyzed the dynamical behavior of branching annihilating ran- 
dom walks in one dimension. Typically this system exhibits (i) propagation, where a 
localized population of particles spreads ballistically for large branching and/or small re- 
action rates, and (ii) extinction, where a localized population eventually disappears for the 
complementary range of parameters. For parity non-conserving rules, there is a continu- 
ous transition between these two behaviors which occurs at a finite value of the branching 
rate when the reaction probability is unity. This transition appears to have a universal 
character and can be accounted for by analysis of a truncated hierarchy of rate equations 
for multiparticle distribution functions, as well as by scaling. 

When the branching and reaction mechanisms conserve parity, rather different be- 
havior arises. In particular, for the 2-offspring BAW with a unit reaction probability, an 
exact analysis of the master equation for the parity of the particle number in a fixed length 
interval shows that the transition between propagation and extinction occurs for infinite 
branching rate. The long-time behavior also depends fundamentally on the parity of the 
initial number of particles. When this number is odd, extinction corresponds to a long 
time state where the particle number is bounded, while the particle number vanishes as 
t- 1 ' 2 for an even number of initial particles, i. e., extinction is complete. In the limit of 
infinite branching rate, a single initial particle gives rise to a Gaussian-like density distri- 
bution where the particle number grows as t 1 ! 2 . Numerical studies confirm many of these 
analytical predictions. 
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Figure Captions 



Figure 1. Phase diagram of the BAW process as a function of the reaction probability k and the 
ratio of the diffusion to branching rate D/r. In (a) the phase diagram appropriate 
for BAW processes with non-parity conserving branching is indicated. This behavior 
should be taken as representative of the continuum limit of the BAW process. In (b), 
the phase diagram for the 2-offspring BAW is schematically indicated. In this case, 
the nature of the "extinct" phase depends on whether the initial number of particles 
is even or odd. 

Figure 2. Illustration of the space-time evolution of the 2-offspring BAW in one dimension in the 
pure branching limit, D/r = 0, and k = 1. In this case, an isolated nearest-neighbor 
pair of particles remains bound and executes a simple random walk. The particle 
which branches at the next step is indicated by the open circle. 

Figure 3. Enumeration of the processes which can change the parity of the number of particles 
contained in the interval [j, k] under the action of the 2-offspring BAW. Shown are 
the various microscopic events and their corresponding statistical weights. The cases 
of (a) k > j and (b) k = j are somewhat different and therefore need to be treated 
separately. 

Figure 4. Enumeration of the processes which can change the occupancy of a selected set of sites 
in the 1-offspring BAW. In (a) the set consists of a single site, po, while in (b) the 
set consists of 2 adjacent sites, po,i- In (c), the set consists of 2 non-adjacent sites, 
Po,i- The signs to the left denote whether a process increases (+) or decreases (— ) p. 
The rate of the processes (and their kind) is indicated to the right. Solid and empty 
circles denote occupied and empty sites, respectively. 

Figure 5. Spatial distribution for the l-offspring BAW in one dimension at t = 1.5 19 = 2217. 

The cases shown are: (a) 25000 configurations at D/r = 1/.88, where the survival 
probability has decreased to approximately 0.11, and (b) 200 configurations at D/r = 
1/.95, where the survival probability has saturated at approximately 0.81. 
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